Resonance mode in B\ g Raman scattering — a way to distinguish between 
spin-fluctuation and phonon-mediated d— wave superconductivity 
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We argue that Raman scattering in B\ g symmetry allows one to distinguish between phonon- 
mediated and magnetically-mediated d— wave superconductivity. In spin mediated superconductors, 
Big Raman intensity develops a resonance at a frequency u> res < 2A max , whose origin is similar 
to a neutron resonance. In phonon-mediated d— wave superconductors, such a resonance does not 
develop. Several extensions of the argument are presented. 



This paper is devoted to the analysis of whether there 
exists an observable that would distinguish between mag- 
netically mediated and phonon-mediated d— wave super- 
conductivity. This subject is motivated by the cuprates 
for which recent measurements, particularly the obser- 
vation of the kink in the quasiparticle dispersion 
has revived the discussion of whether the pairing in the 
cuprates may be due to phonons rather than spin fluc- 
tuations. Both clectron-phonon |2| and spin-fermion in- 
teractions were advanced to explain the features in 
the quasiparticle dispersion. Thus, to truly distinguish 
between the two scenarios, one needs to identify an ob- 
servable for which phonon and spin-fluctuation pairings 
yield qualitatively different results. 

In this communication, we argue that Raman scat- 
tering in Big symmetry is such a probe. We show 
that in spin mediated d— wave superconductors, the B\ g 
Raman intensity develops a resonance at a frequency 
while in phonon-mediated d— wave su- 
perconductors, such a resonance does not develop. The 
resonance is excitonic in origin, and is similar to the ex- 
citonic resonance in the spin susceptibility of a d— wave 
superconductor p|. The only difference is that the reso- 
nance in B\g Raman scattering comes from fermions all 
around the Fermi surface and has a finite intrinsic width 
because of the nodes, whereas the resonance term in the 
spin susceptibility comes only from fermions in the antin- 
odal regions and is a true, 8— functional bound state. 

To make the argumentation straightforward, we first 
consider the simplest case of S = 1/2 fermions interacting 
via a static potential V patr (k) 

Hint = - 53 ^L^p+q,/3 VZi&sQt-P) ^P7^k+q,5. (1) 

q,k,p 

In Eq.(l) summation over spin indices a, [3, 7 and 8 is 
understood. As B\ g vertex has the same d— wave form 
as the pairing gap, we further approximate V pair by its 
d— wave component V pair {k — p) oc dkd p , where dk = 
[cos(k x a) — cos(k y a)}/2. 

The effective interaction V p ^ s (k — p) may be due 
to spin fluctuations or to phonons. For spin mediated 



interaction, V p ^ s (k - p) = V spin d k d p (J Q7 ■ crps where 
a are Pauli matrices. For phonon-mediated interaction 
V^g*I 5 (k— p) = Vphdkd p S ai S/3s- We assume that both in- 
teractions lead to a d— wave pairing, and study the conse- 
quences for the Raman intensity. ^From this perspective, 
our results are equally applicable if phonons are replaced 
by charge density waves [5j . 

Our first observation is that the signs of V sp i n and 
Vcharge must be different, if they both lead to an at- 
traction in a d— wave channel. Indeed, substituting effec- 
tive interctions into the diagrammatic expression for the 
d— wave, spin-singlet pairing vertex tp^ a&difyapTpLu r> 
where $d{k) a p = ^d k a v al3 , and using (T v a pS ai 8ps = 
a y8> (jV afi' J a- 1 ■ crps = -3er^, we obtain that <& is re- 
lated to the bare vertex <3>o as 

d>= ?2 ; cj> = ( 2 ) 

1 + 2>AV spm 1 - AV p h 

for spin-mediate or phonon-mediated interaction, re- 
spectfully. A positive A cx | log ui c \ is a conventional log- 
arithmic cut-off factor. To get an attraction, one then 
obviously needs V p h to be positive, V sp in to be negative. 
This is the case when phonon mediated interaction is 
peaked at small momenta, and spin-mediated nteraction 
is peaked at momenta near (ir,pi) 0,0- 

Suppose next that the system is a d— wave supercon- 
ductor, either due to phonons or due to spin fluctua- 
tions. Consider how this affects the Raman response in 
the Big channel. The Bi g Raman vertex T a p(k) has 
d— wave momentum dependence and is a spin scalar: 
r a p(k) = TdkSap. For a BCS superconductor without 
vertex or self energy corrections, the Bi g Raman intensity 
is lg (D) = r 2 7mxo(^) and Xo(^) is the imaginary part 
of the particle-hole bubble with two d-wave vertices (see 
e.g. Ref.Q). In the normal state, Xo(^) vanishes, as 
it should as there is no low-energy phase-space available 
for scattering with q = 0. In the superconducting state, 
light can break Cooper pairs with q = and Yn(f2) is 
given by Tsuneto function weighted with d\ [8|, |s| The 
imaginary part of Xo(^) scales as Q 3 at small frequen- 
cies [9|, and diverges logarithmically as f2 approaches 
±2A: Imxo(^) oc log[A/V^ 2 - 4A 2 ]. 
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FIG. 1: The real part of the B\ g Raman response in the 
superconducting state Rexo for the case of small fermionic 
damping (solid line) and no fermionic damping (dashed line) . 
The vertical dotted line marks 2Ao. A small change between 
2Ao and the frequency where Rexo(Q) is peaked is due to 
the fact that the maximum value of the gap along the Fermi 
surface is slighly smaller than Aq. 



The corresponding real part at small frequencies varies 
as Rexo(^) = N [l + 0(tt 2 /A 2 )], with N the density of 
states at the Fermi level, and keeps increasing up to 2A 
before discontinuously jumping across zero at a frequency 
of twice the maximal gap on the Fermi surface. An- 
other milder jump occurs at twice the energy at the van 
Hove points at (jr, 0) and related points in the Brillouin 
zone. We plot Re%o(^) m FigQ] using A = 35 meV and 
the band structure e k given by Eschrig and Norman [To) 
to fit the angle-resolved photoemission (ARPES) on op- 
timally doped Bi 2 Sr 2 CaCu 2 8+2; (Bi-2212). Although 
the generic behavior changes near 2Ao if damping (phe- 
nomenologically represented by 8) is increased, the rapid 
rise and fall of the real part near 2A is preserved. 

Consider now what happens to the Raman response 
when we add V pa i r . This interaction leads to two ef- 
fects - self-energy renormalization of the fermions in the 
particle-hole bubble, and the renormalization of the Ra- 
man vertex. The self-energy renormalization does not 
distinguish qualitatively between phonons and spin fluc- 
tuations, and still preserves the peak Ib 1s (ft) at ft = 2A. 
The renormalization of the B\ g vertex is more relevant. 
To understand its role we first observe that there is no 
spin-induced sign change between vertex renormaliza- 
tion due to phonons and due to spin fluctuations. In- 
deed, convoluting the spin dependence of V^ l ^ s with 
Sp-y of the Raman vertex T({3"f)(k), we find that for 
magnetic interaction, the summation over spin indices 
yields b s ^a ai ap S = J2i( a as)i = 3$aS, while for phonons 
Sfj-fda-ySfjs = S a s- As a result, in both cases, the spin 
configuration of the Raman vertex is reproduced, and 
the sign is the same in both cases. Summing up the ver- 
tex correction diagrams in the ladder approximation, we 



obtain 

jph /^\ _ p2 Xoffi) 

Bla (i + iv spmX ' (n)y + (lv spmX '^n)Y 

We defined x' — R e X, x" — Imx- We now recall that 
V S pin must be negative, and V p h must be positive, if each 
of them gives rise to a d— wave pairing. Then, the sign of 
the vertex renormalization in is different for phonons 
and spin fluctuations. As mentioned above, Im Xoity m 
the superconducting state is quite small except for near 
2 A. The renormalization of the Raman vertex at fl < 2 A 
then predominantly comes from Kexo- 

Since Rexo(^) is positive, vertex renormalization due 
to phonons reduces the Raman vertex at small frequen- 
cies and only slightly shifts up the peak which remains 
close to 2A. On the other hand, if the d— wave inter- 
action is magnetic in origin, V sp i nX ' (fl) < 0, and for 
strong enough V sp i TL , there exists a frequency Q re s < 2 A 
at which (3/2)V sp i n x'o(Q) = —1- At this frequency, 
jspm a pggjj. Because Xo(^) i s nonzero at any 
O > 0, the peak is not infinitely sharp as in an s— wave 
superconductor^. However, since Xo ^ 3 at small 
frequencies, the width of the peak is small if fl re s is sub- 
stantially smaller than 2A. We see therefore that for 
spin-mediated d— wave pairing, the B\ g Raman intensity 
develops a resonance. The resonance for spin-mediated 
pairing was discovered in Ref , although its origin was 
not discussed in detail there. We plot in Figure|21the full 
Big Raman response for T = in the superconducting 
state for different values of the interaction N V using the 
parameters shown for Figure In the absence of inter- 
actions (V = 0) the Raman response rises as f2 3 and has 
a clear peak at twice the gap and another smaller peak 
at twice the van Hove energy. For magnetic interactions 
V < the low energy peak sharpens and moves to lower 
frequency as the resonance develops and steals spectral 
weig ht from the 2Ao feature. In agreement with Ref. 

two separate peaks do not develop, but because 
of resonance, the original 2A peak progressively shifts 
down from 2A as the interaction increases. Conversely, 
for phononic interactions V > 0, the peak renormalizes 
upwards from the final-state interactions and weakens for 
stronger interactions. No low energy resonance develops. 

The resonance in the B\ g Raman response is in many 
respects similar to the resonance in the spin susceptibil- 
ity in a d— wave superconductor Q. In both cases, the 
d— wave symmetry of the gap is crucial, and the resonance 
emerges due to residual attraction between fermions and 
spin-fluctuations in a d— wave superconductor. The only 
difference between the two resonances is that the neu- 
tron resonance is virtually a bound state in the sense 
that it is infinitesimally narrow, and does not require 
coupling to be above a threshold because the real part of 
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FIG. 2: The full Raman response plotted for VN = 
(black solid line), and ±0.01, 0.1, 0.5 (dotted, dashed, dashed- 
dotted lines, respectively) for phonon-mediated (green lines) 
and spin-mediated (red), respectively. 



the bare spin susceptibility Xs{Q,&), where Q = {n,ir), 
evolves between and infinity at < £1 < 2 A, where 
Im Xs(Qity — 0. This comes about because the spin 
susceptibility at momentum Q is determined by fermions 
near hot spots, where k and k + Q are both near the 
Fermi surface. The hot spots are generally located away 
from the nodes, hence in the superconducting state hot 
fermions are fully gaped, and a spin fluctuation needs 
a finite energy to be able to decay into a particle-hole 
pair. The Raman resonance, on the other hand, is a 
Q = probe. It therefore involves fermions from the 
entire Fermi surface, including nodal regions. The nodal 
fermions account for a nonzero Imxo(^) at any finite 
and therefore produce an intrinsic width of the Raman 
resonance peak. The interaction then should be above 
the threshold for the resonance to become visible. 

We now extend the above analysis in several different 
directions. First, when V sp i n and V„h are both nonzero, 



ph 



the resonance condition becomes (l/2)V e //Xo(^) = — 1; 

V p h- We remind the reader that 



with V e ff = 3Vs 



spin 



Vspin < and V p h > both favor d— pairing, and thus 
they compete to determine V e ff. Thus referring back 
to Figure (0) the full B\ g Raman response is a function 
of the net pairing interactions V e ff, and the presence of 
the resonance requires 3V sp in >| Vphonon | if d— pairing 
occurs via both channels. 

Second, the analysis can be extended to the case where 
the pairing interaction and the Raman vertex have the 
same symmetry, but are not identical. If V pair (k — p) oc 
'/'k^p and the bare Raman vertex is Tk = T7k, where the 
functions (/>k and 7k belong to the same irreducible rep- 
resentation different from fully symmetric A\ g , the RPA 
Raman intensity near the resonance at a temperature T 
is given by 



7^(0) ~ Im 



& 2 K w K m 



1 + 2A 2 (3K 
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V ph )K { 
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iSf) tanh 2T- 
(5) 



One can verify that for 



d k , K 00 (Cl,0) = #io(0,0) = #oi(fi,0) = (l/4Ag) X o(fi), 
and (@J coincides with Eq. (@J however shows that 
once SVspin + V p h is negative, the resonance occurs for 
any pairing symmetry different from A\ g , for the Raman 
vertex for which Kqi and Kiq are finite. 

Third, we argue that for spin mediated pairing d— wave 
pairing, the resonance in the B\ g channel survives even 
when the effective interaction includes a dynamic term, 
e.g., the Landau damping, and vanishes at high frequen- 
cies. The vanishing of V pair (k, uj) at large ui may have 
a profound effect on the renormalization of the vertex 
at zero momentum transfer as the real part of the ver- 
tex renormalization, which accounts for the resonance in 
our case, partly comes from fermions with high frequen- 
cies. Once the interaction vanishes at high frequencies, 
this part disappears, and it becomes an issue whether the 
remaining part still has the same sign. 

To address this issue, we computed V sp i n xo(ty i n 
Eliashberg theory which includes Landau damping. We 
obtained the expression for the Raman vertex renormal- 
ization rs ls (S!) = r<ifc/(l + J(f2)) by first evaluating 
the vertex correction in Matsubara frequencies, and then 
analytically continuing to the real axis by introducing 
double spectral representation. We found that J(0) = 0, 
but J(f2) < at f2 < 2A still logarithmically diverges at 
2 A, i.e., vertex corrections, even without high frequency 
term, still increase the Raman vertex and lead to a reso- 
nance in R(£l) at some ft < 2A. 

The only extensive set of data comparing B\ g sym- 
metry Raman gap values with those from the single 
electron spectroscopies of angle-resolved photoemission 
spectroscopy (ARPES) and tunneling is found for the 
Bi 2 Sr 2 CaCu20 8+x ("Bi-2212") family. Values found by 
five groups for the Raman gap in terms of the hole dop- 
ing p are shown in Fig [3] The Raman gap values are 
compared in this figure with those from tunneling and 
ARPES. Tunneling results are from Ref [24| and repre- 
sent peak-to-peak separations in positive and negative 
biases. ARPES results are from Ref For complete- 
ness, we presented twice the gap A determined from two 
sets of ARPES data - the position of the peak at (0,7r), 
and the mid-point of the leading edge gap inferring from 
several different forms of modelling. 

For doping p greater than 0.2, well-defined peaks 
emerge below T c in the B\ g channel at a frequency 
roughly consistent with the tunneling gap. Both Raman 
and tunneling data fall below rather scattered ARPES 
data The agreement between Raman and tunneling re- 
sults indicates that the B\ g vertex renormalization is 
small in the overdoped regime. This may be the con- 
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FIG. 3: The Raman data (open squares) from five different 
groups (Ref. Il3| (red), Ref. |2Cl (magenta), Ref. |2l| (or- 
ange), Ref. |3T(green), Ref. [l| (blue), Ref. [l<| (cyan), 
respectively) , together with the tunneling data (red triangles) 
from and ARPES data from [g3 (black circles are de- 
rived from leading edge spectra while green circles are from 
the peak position of the 7r,0 spectra). The doping p was de- 

termined by the formula p — 0.16±i / ^ 6T c ma:c , and the solid 

line plots 70 meV xT c /T™ ax . 



from Big Raman, ARPES, and STM and other probes 
would not be unexpected. This remains very much an 
open question meriting further investigation. 

To conclude, in this paper, we considered B\ g Ra- 
man intensity in a d— wave superconductor. For non- 
interacting fermions, B lg intensity is peaked at 2A. We 
found that for interacting fermions, there is a qualitative 
distinction between spin-mediated and phonon-mediated 
interactions. For spin-mediated interaction, the peak in 
the intensity shifts downwards compared to 2A due to 
the development of the resonance below 2 A. We verified 
that the resonance survives even if the interaction is re- 
tarded. For phonon-mediated interaction, the resonance 
does not develop, and the peak remains at 2A. The effect 
is quite generic and is also valid if phonons are replaced 
by charge-density-waves. We presented the generic con- 
dition under which resonance occurs in a Raman vertex 
of arbitrary symmetry. 
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and D. Morr for useful discussions. The research was 
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von Humboldt Foundation, ONR Grant No. N00014-05- 
1-0127, and NSERC (T.P.D.). 



sequence of just small enough spin-fermion coupling, or 
the partial cancellation between the renormalizations due 
to spin fluctuations and due to phonons. The distinction 
with ARPES is likely due to ARPES resolution and also, 
possibly, to the effects from bilayer splitting. 

For p less than or equal to the optimal value, the Ra- 
man gap values in Fig [3] are mostly from Ref i 
(l9|. and Ref |2jj. They fall consistently below those 
from tunneling and from ARPES although there is some 
degree of scatter. In regard to the findings of this paper, 
this implies that the pairing more likely arises from spin 
fluctuations in this range of doping. Yet the growth of 
intensity predicted in Figure @ is not seen. This may be 
a result of inelastic scattering or may be due to the pres- 
ence of the pseudogap. In addition for reasons that are 
not well understood, other groups 0, |2^ find very weak 
Big Raman signals from underdoped samples of Bi2212. 
They thus find it difficult to impossible to extract values 
of 2 A in this symmetry and doping range. 

Yet this is not the case for £?2 S scattering (light orien- 
tations rotated by 45 degrees with respect to Bi g orien- 
tations). Peaks have been shown to occur only at tem- 
peratures below T c at a frequency with scales with T c 
for all dovinas\ll\ . This apparent contradiction between 
the findings in Bi g and Bi g requires more study and 
must be reconciled before firm conclusions can be drawn 
about the pairing mechanism. A possible alternative sce- 
nario could be that the anti-nodal quasiparticles become 
gapped via a mechanism which is not related to super- 
conductivity, such as precursor spin- or charge-density 
wave pairing. In this case, the difference in values of 2A 
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